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摘 要 : 动态 自 适应 多 尺度 小 波 配点 法 (AWCM) 能 有 效 地 模拟 具有 间歇 性 的 物理 现象 ， 此 方法 是 近 几 年 发 展 起 来 的 非常 新 颖 


的 数值 计算 方法 。 为 了 增强 该 方法 识别 与 跟踪 解 的 奇异 性 的 能 力 ， 并 提高 数值 计算 稳定 性 与 计算 效率 ， 本 文 将 精细 时 程 积分 算 
法 与 之 相 结 合 形成 了 快速 动态 自 适 应 多 尺度 小 波 配 点 法 。 为 了 实现 这 一 算法 , 本 文 给 出 了 构造 动态 自 适应 网 格 配 点 集 的 新 方法 ， 
构建 了 以 小 波 ( 或 尺度 函数 〉 系数 为 变量 的 时 程 推进 公式 。 通 过 求解 一 维 Burgers 方程 ， 证 明了 本 文 方法 具有 更 加 良好 的 数值 
计算 性 质 。 
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A Fast Dynamically Adaptive Wavelet Collocation Method 


Based on the Precise Time-integration 
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Abstract: The dynamically adaptive wavelet collocation method (AWCM) is a new method developed recently for 
simulating the physical phenomenon with intermittency efficiently. In order to enhance the stability and efficiency of 
computation and the ability of identifying and tracking singularity, a new method is established combined the precise 
time-integration method with AWCM. In this method, a new way to generate dynamic adaptive grids is proposed and 
the time marching equations as functions of wavelet (or scale function) coefficients are set up. Validated by solving 
one- dimensional Burgers equation, it ls found that the new method provided in this paper presents good 
performance of computation. 
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0 前 言 

在 计算 流体 力学 方面 ， 小 波 方法 构成 了 一 个 非 动态 自 适应 多 尺度 小 波 配 点 法 四。2000 年 ， 在 第 二 
常 新 颖 的 研究 领域 。 此 方法 最 早 可 追溯 到 1990 年 ， 代 小 波 的 基础 上 构建 了 求解 偏 微分 方程 4 第 二 代 小 
Glowiski 等 由 使 用 小 波 方法 成 功 求解 椭圆 型 与 双 曲 波 配点 法 四。 在 小 波 方 法 求解 满 流 流动 方面 ,最 早 使 
型 偏 微分 方程 。Vasilyev 在 此 方面 做 出 了 巨大 的 页 ” 用 的 方法 是 基于 小 波 的 直接 数值 模拟 (WDNS)， 此 
献 。1995 年 ， 在 求解 偏 微分 方程 配点 法 与 小 波多 分 。 方法 充分 考虑 了 满 流 流动 的 间 欣 性 来 降低 计算 复杂 
辨 分 析 的 基础 上 , 提出 了 多 尺度 小 波 配点 法 中。1996 度 。 其 次 是 1999 年 Farge 等 中 提出 的 相干 涡 模 拟 方 
年 ， 在 其 提出 的 多 尺度 小 波 配 点 法 的 基础 上 深入 看 法 〈CVS)， 此 方法 利用 小 波 将 淇 流 流 动 分 解 为 拟 序 
究 ， 利 用 小 波 所 具有 的 信号 奇异 性 检测 能 力 构 建 了 
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部 分 与 随机 背景 流动 部 分 ， 前 者 直接 计算 ， 而 后 者 
建立 简化 模型 求解 。2004 年 ，Goldstein 等 四 提出 了 
随机 相干 涡 自 适应 大 涡 模 拟 〈SCALES)， 此 方法 继 
承 了 相干 涡 模 拟 的 跟踪 庙 流 流动 高 含 能 拟 序 结构 的 
能 力 ， 同 时 使 用 了 大 涡 模 拟 (LES ) 模 化 低 含 能 随机 
i 
1994 年 ， 钟 万 砚 思 院士 提出 了 结构 动力 学 的 精 
时 程 积分 法 ， 此 算法 只 在 空间 上 使 用 近似 的 离散 
格式 ， 对 于 线性 问题 此 方法 是 一 种 半 解 析 数 值 模拟 
方法 ， 通 过 此 算法 对 时 层 推进 使 用 的 常 指 数 和 矩阵 进 
行 精细 计算 ， 使 其 在 时 间 上 具有 很 高 精度 
本 文 综合 考虑 数值 计算 效率 与 稳定 性 ， 提 出 了 
种 基于 精细 时 程 积 分 的 快速 动态 自 适 应 多 尺度 小 
波 配点 法 。 
1 计算 方法 
本 文 所 采用 的 是 基于 插值 小 波 变换 的 动态 自 适 
应 多 尺度 小 波 配点 法 。 由 于 小 波 变换 只 能 在 能 量 有 
限 的 函数 空间 ZJ 进行 ， 因 此 本 文 只 针对 ZJ 空间 
函数 展开 讨论 。 以 一 维 Burgers 波动 为 例 进行 介绍 
其 无 量 纲 控 制 方程 为 : 
Gu Ou_ 1 O%u 
ot a Re Ox 
1.1 插值 小 波 变换 离散 偏 微分 方程 
对 u(x) 进 行 如 下 小 波 变 换 : 
2 


4) Su bn Ee ld) 


如 =0 />j0 k=0 


AS 


(1) 


原 函数 u(x) 关于 x 的 各 阶 导数 可 以 由 尺度 函数 
p(k) 与 小 波 函数 y (k) 的 各 阶 导数 来 计算 ; 
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uy (x) = Dl ul jh ja + py oy i (x) . 


ko=0 jjo k=0 (3) 


将 式 (2)、 式 (3) 代 入 方程 (1) 得 到 一 维 Burgers 
的 半 离 散 方 程 : 


Ou J (xn,t) 


2 如 
| 
ot 隐 
2 7]0 


we Dwr bn (1)+ > So (%) 


=0 六 jn k=0 


Xi | Xn) + Sk 


Jj2jo k=0 


(4) 
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最 细 尺 度 ; 


点 上 的 近似 值 ; 
u(xin ) 一 一 最 粗 尺度 配点 集 的 第 个 配置 点 


尺度 函数 Gi (x) 的 系数 ; 


pit (x ) i 集 的 第 如 个 配置 点 


尺度 函数 p，(x) 在 最 细 尺 度 配点 集 的 第 个 配置 点 


最 粗 尺度 配点 集 的 第 个 配置 点 


尺度 函数 办 ，(x) 在 最 细 尺 度 配点 集 的 第 个 配置 点 


上 二 阶 导 
9 第 大 个 配置 点 小 波 函数 
yi (zx) 的 系数 ; 
ya) 集 的 第 个 配置 点 小 波 


函数 yi (x) 在 最 细 尺 度 了 配点 集 的 第 n 个 配置 点 上 
一 阶 导 数 ; 


Wn (mw ) 


集 的 第 上 个 配置 点 小 波 


函数 yi (x) 在 最 细 尺 度 配 点 集 的 第 n 个 配置 点 上 


二 阶 导数 。 
1.2 基于 精细 时 程 积分 的 多 尺度 插值 小 波 配 点 法 
半 离 散 方程 《4) 可 以 写成 如 下 矩阵 形式 : 
SU bg C—diag(U)SC (5S) 


ot Re 
式 中 ， 向 量 C 由 小 波 ( 或 尺度 函数 〉 系数 构成 ; 向 
二 阶 导 


S; 分 别 由 小 波 〈 或 尺度 函数 ) 一 阶 、 


量 $i、 
数 构 成 。 
半 离 散 方程 (5) 中 向 量 T 与 向 量 C 存在 下 式 关系 : 

U=S0C (6) 


向 量 So 由 小 波 《〈 或 尺度 函数 ) 构成 。 


矩阵 So 是 满 秩 矩 阵 ， 存 在 逆 和 矩阵 So ,可 得 : 


C=SJU (7) 
将 式 (7) 代入 式 (5) 可 得 : 
全 -HU+diag(U)MU (8) 


式 中 ，H -了 SS ，M =SISi 。 
Re 


方程 (8) 构成 了 U 关于 上 的 常 微分 方程 ， 具 有 
解析 解 : 


U(t, ) 三 exp| Hu 一 塘 1 )|U(#) 


-exp(Hi)| diag[u(O]Mu()d 


(9) 


等 式 (9) 右 端 第 二 项 是 Burgers 方程 的 非 线 性 项 ， 
为 了 实现 积分 需要 进行 线性 化 处 理 ,在 t= 处 Taylor 


展开 : 


diag| U(7) [MU()=rotn (t -ti) 
m=diag| U(t) |MU(i,) 
se wo (10) 
Ot | 
ee 
OU 
diag| U(t,) IM| 一 一 
十 iag| (i)| 攻 | 


将 式 (10) 代入 式 (9) 并 整理 可 得 方程 (8) 
的 时 层 推进 公式 : 
U(s)=TU(i)+(T-DH (n+H) ee 
-Hi 人 si fi) 


其 中 ， 指 数 矩 阵 T= exp[HC -1)] 可 以 精确 求解 
首先 将 7 表示 成 以 下 形式 : 


T= exp|HC 对 ti 1)] = exp[H ry 
四 At ft (12) 
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取 M =20，z 将 是 一 个 非常 小 的 时 间 量 ，Hzr 接 


近 于 零 矩 阵 。 将 exp[Hz] 在 零 矩 阵 处 Taylor 展开 : 


exp[H7| =]+Ta 


=I+Hr+(Hr) [I+(H7)/3+(H7) /2|f2 


(13) 


值得 注意 的 是 Ta 接近 零 第 阵 ， 必 须 与 单位 阵 I 
分 开 存 储 以 保证 其 精度 ,将 Ta 按照 下 式 进 行 M 次 迭 
代 ; 

Ta, =(I+Tai) -1=2Ta +Ta? 


Tas =(I+Ta) -1=2Ta, +Ta3 


Tav =(I+Tay 1) -7=2Tay i+Tay ， 


可 以 得 到 T 和 矩阵 的 高 精度 解 T=I+ Tayv 。 


1. 3 时 层 推进 的 动态 自 适 应 网 格 的 新 方法 
在 传统 的 动态 自 适应 多 尺度 小 波 配 点 法 实施 过 
程 中 ， 为 了 使 数值 计算 能 够 跟踪 解 的 奇异 性 ， 不 仅 


保留 了 已 求解 的 上 一 时 层 上 满足 条 件 


j 1/2 
Cl2al é 


《a2 为 与 尺度 有 关 常 数 ， 为 小 波 系数 阔 值 ) 的 配 


置 点 ， 还 要 保留 待 求解 的 当前 时 层 可 能 对 解 的 奇异 
性 有 贡献 的 重要 配置 点 。 由 于 数值 计算 的 时 间 步 长 
较 小 ， 解 的 奇异 性 只 能 从 上 一 时 层 配 置 点 迁移 到 其 
附近 配置 点 ， 因 此 ， 当 前 时 层 上 可 能 对 解 的 奇异 性 
有 贡献 的 重要 配置 点 位 于 上 一 时 层 满 足 条 件 


cj|> ae (af 为 与 尺度 有 关 常 数 ，s 为 系数 阔 值 ) 


的 相 邻 配置 点 。 
对 传统 的 动态 自 适应 多 尺度 小 波 配 点 法 的 不 规 
则 配点 集 的 构造 方法 进行 分 析 ， 并 通过 数值 试验 进 
行 验证 ， 给 出 以 下 结论 : 当时 间 步 长 At 过 大 ， 传 统 
的 动态 自 适 应 多 尺度 小 波 配 点 法 不 能 跟踪 解 的 奇异 
性 ， 其 不 规则 配点 集 的 构造 方法 是 影响 数值 计算 稳 
定性 的 重要 原因 。 

为 了 增强 数值 计算 的 稳定 性 ， 本 文 改进 了 动态 
自 适 应 网 格 构建 方法 一 一 边 求解 边 生 成 的 方法 。 在 基 
于 精细 时 程 积分 的 快速 动态 自 适应 多 尺度 小 波 配点 
法 的 实施 步骤 中 给 出 此 方法 : 
(1) 选择 插值 小 波 ， 构 造 时 层 推进 所 需 的 矩阵 A、B 
与 Si; 


(2) 假设 在 不 规则 配点 集 G2 上 获得 了 上 一 时 层 的 


尺度 函数 或 小 波 函 数 系数 c 站 ， 计 算 1 丰 时 层 w 函数 


在 配点 集 G 上 的 近似 解 ur; 


(3) 在 配点 集 G 中 最 粗 尺度 的 尺度 空间 配点 集 G。 


上 ， 计 算 当前 时 层 系 数 cec ; 
(4) 假设 配点 集 G 中 较 小 尺度 的 小 波 空 间 配 点 集 
G, 与 G,、G 合 并 形成 G+ ， 在 G; 中 G 上 的 配置 点 


与 G,、G 中 的 配置 点 相间 构成 临 点 。 若 G, 中 配置 


点 相 邻 配置 点 的 系数 


ceaua|> ajps , 则 计算 此 配置 


点 上 的 系数 cee 


(5) 与 步骤 〈4) 相 类 似 的 计算 其 他 小 波 空间 配点 
集 上 的 系数 cv ， 最 终 在 配点 集 G 上 获得 当前 时 层 可 


能 对 解 的 奇异 性 有 贡献 的 配点 上 的 系数 c, ; 


(6) 通过 式 |c|> ay2e 保留 配置 点 ， 构 造 当前 时 层 


的 不 规则 配点 集 G4 ， 计 算 上 时 层 z 函数 在 配点 集 G 
上 的 近似 解 2 
(7) 重复 (2) ~ 〈6) 进行 时 层 推进 
终 时 层 1 上 函数 的 近似 解 w*”。 
2 算 例 分 析 与 讨论 
2.1 问题 描述 

方程 (1) 给 出 了 一 维 无 量 纲 Burgers 方程 ， 计 
算 区 域 为 (0,2) ，Re=1000， 给 定 其 初始 条 件 与 边界 


直至 


CX. 
: 
川 


条 件 为 : u(x,0) = sin(xx)， 


采用 在 时 间 上 具有 一 阶 精度 的 显 式 格式 的 动态 
自 适应 多 尺度 小 波 配 点 法 与 一 阶 精度 的 基于 精细 时 


u(0,0D) =u(2,1)=0。 
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性 的 物理 量 场 ， 局 部 的 奇异 点 意味 着 局 部 的 高 频 波 


动 ， 导 致 细 尺度 配置 点 上 小 波 系 数 增 大 ; 随 着 数值 
解 局 部 奇异 性 的 增强 ， 数 值 计算 能 够 相应 增加 细 尺 
度 配置 点 ， 体 现 了 算法 具有 识别 奇异 点 的 能 力 。 
1.0 220 7F 
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sl 
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图 1 Burgers 方程 的 解 与 配置 点 的 时 间 演 化 过 程 
Fig.1Time evolution of solutions of Burgers equation and 
collocation points 

为 了 对 比 两 种 算法 的 稳定 性 ， 在 时 间 步 长 
A 太 10”、A 广 10 ， 小波 阔 值 =10” 条件 下 进行 数值 模 
拟 。 图 2 给 出 了 一 维 Burgers 方程 使 用 两 种 算法 在 各 
时 间 步 长 下 的 计算 结果 与 配置 点 使 用 情况 。 从 图 2 
中 可 以 看 出 ， 在 A 记 10*、Af=10” 时 间 步 长 下 ， 本 文 
方法 都 能 够 有 效 地 识别 、 捕 捉 与 跟踪 解 的 奇异 点 ， 
并 对 物理 问题 进行 精确 模拟 ， 传 统 方法 在 较 小 时 间 


程 积 分 的 快速 动态 自 适应 多 尺度 小 波 配点 法 进行 求 
解 ; 边界 处 采用 对 称 延 拓 ; 将 计算 区 域 均匀 离散 为 
2+1 个 配置 点 。 
2. 2 计算 结果 分 析 

设 定 小 波 阔 值 =10”， 信 t=10-4， 使 用 本 文 提出 
的 基于 精细 时 程 积分 的 快速 动态 自 适 应 多 尺度 小 波 
配点 法 对 上 述 两 类 物理 问题 进行 求解 ， 图 1 给 出 了 
各 无 量 纲 时 间 上 的 数值 计算 结果 ， 以 及 数值 计算 所 
使 用 的 动态 网 格 配置 点 。 观 察 图 1，Burgers 方程 从 
光滑 的 初始 物理 量 场 逐渐 发 展 为 局 部 具有 强烈 奇异 


A 


步 长 A 万 10 条件 下 ， 同 样 可 以 精确 模拟 物理 过 程 ， 
然而 在 较 大 时 间 步 长 A 二 107 条 件 下 ， 随 着 解 的 奇异 
性 增强 ， 误 差 逐渐 积累 放大 ， 当 大 0.38 时 ， 开 始 产 
生 振 荡 ， 数 值 计算 发 散 。 上 述 结果 充分 说 明了 本 文 
提出 的 新 算法 具有 更 高 的 数值 计算 稳定 性 ， 其 原因 
为 : 尽管 两 种 算法 在 时 间 上 都 具有 一 阶 精度 ， 但 
是 本 算法 在 时 层 推 进 方面 是 基于 精细 时 程 积分 方 
法 ， 能 有 效 抑制 误差 的 放大 ;@ 传 统 方法 不 规则 配 
点 集 的 构造 方法 必须 要 求 较 小 的 时 间 步 长 ， 否 则 将 
失去 跟踪 解 的 奇异 点 的 能 力 ， 而 本 文 算法 的 不 规则 


配点 的 构造 对 时 间 步 
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(al) A 六 10“， 本 文 算法 


(a2) A 太 103， 本 文 算法 
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2 两 种 算法 计算 结果 与 配置 点 


Fig.2 Comparisons of Solutions and collocation points using two 


method 
表 1 给 出 了 两 种 算法 在 各 时 间 步 长 与 小 波 阔 值 
条 件 下 CPU 时 间 ， 从 表 1 中 可 以 看 出 ， 传 统 算法 在 
较 高 时 间 步 长 下 计算 发 散 ， 在 较 低 时 间 步 长 下 CPU 


时 间 远 高 于 本 文 方法 ， 其 计算 时 间 是 本 文 算法 的 10 
倍 左 右 。 
表 1 计算 时 间 对 比 
Table 1 Comparisons of CPU time 
、 ， 本 文 算 法 传统 算法 
计算 时 间 5 可 区 二 
e=10 e=10 E=10 E=10 

和 A 太 10 46.13s 37.83s 420.63s 383.84s 
A 太 10” 5.63s 4.55s 发 散 发 散 
3 结论 


本 文 结合 传统 的 动态 自 适应 多 尺度 配点 法 与 精 
细 时 程 积 分 算法 的 优点 ， 提 出 了 一 个 新 算法 。 为 了 
实现 这 一 算法 ， 本 文 进行 了 如 下 工作 : 
(1) 给 出 了 动态 自 适应 网 格 不 规则 配点 集 构造 
的 新 方法 ， 从 而 使 配点 集 的 构造 不 约束 时 间 步 长 的 
大 小 ， 提 高 计算 稳定 性 ; 
(2) 构建 了 以 尺度 函数 或 小 波 函 数 系 数 为 变量 
的 时 程 推进 公式 。 
数值 试验 证 明了 本 文 方法 具有 识别 、 捕 捉 与 跟 
踪 解 的 奇异 性 的 能 力 ; 通过 数值 计算 对 比 证 明了 本 
文 提 出 的 新 算法 优 于 传统 算法 ， 主 要 体现 在 : 本 文 
提出 的 算法 数值 计算 的 稳定 性 更 强 、 计 算 效 率 更 高 ， 
在 相同 的 条 件 下 ， 计 算 效 率 提 高 了 约 10 倍 。 
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